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Abstract 

We give a complete classification of all static, spherically symmetric solu- 
tions of the SU (2) Einstein- Yang-Mills theory with a positive cosmolog- 
ical constant. Our classification proceeds in two steps. We first extend 
solutions of the radial field equations to their maximal interval of exis- 
tence. In a second step we determine the Carter-Penrose diagrams of all 
4-dimensional space-times constructible from such radial pieces. Based on 
numerical studies we sketch a complete phase space picture of all solutions 
with a regular origin. 



1 Introduction 



In contrast to the coupling of Abelian gauge fields to gravity, the interaction of 
non-Abelian gauge fields with gravity yields a rich structure of interesting solu- 
tions even in the spherically symmetric sector. While the Reissner-Nordstr0m 
black holes arc the only static, spherically symmetric, asymptotically flat so- 
lutions of the Einstein-Maxwell equations, the non-Abelian SU(2) Einstcin- 
Yang-Mills (EYM) theory has a much larger class of such solutions. There 
is an infinite sequence (indexed by the number n of zeros of the Yang-Mills 
potential) of regular 'soliton-type' solutions, known as the Bartnik-McKinnon 
(BK) solutions [1, 2, 3]. Besides there are corresponding families of black holes 
[4, 5, 6, 7, 3] and 'limiting' solutions, when n tends to infinity [8]. 

If a cosmological constant. A, is added to the theory, the asymptotic be- 
haviour of the regular solutions changes from Minkowskian to deSitter for A > 
or Anti-deSitter for A < 0. Furthermore for A > the regular solutions have a 
cosmological horizon. 

Numerical investigations of the static, spherically symmetric EYM equations 
with a cosmological constant have led to a rather clear picture of the most 
important properties of the solutions [9, 10, 11, 12]. Let us briefly recall here 
some of the most relevant numerical results of Refs. [9, 10, 11]. In Ref. [9] 
families of solutions with a regular origin and a horizon with n = 1, 2, . . . zeros 
were found. For < A < Ac(n) the solutions are asymptotically deSitter. At 
A = Ac(n) the global structure of the spatial sections of the solutions changes. 
The area of the 2-spheres has a maximum ('equator') behind the horizon and 
then decreases back to zero, where the solutions have a singularity, implying 
that the spatial sections have the topology of 3-spheres. As A increases the 
horizon and the equator exchange their positions and as A approaches some 
maximal value Ar{n) the radius rh of the horizon shrinks to zero. The solutions 
for A = Ar{n) are globally regular with compact spatial sections of 3-sphere 
topology. The numerical results of Ref. [10, 11] for A < indicate the existence 
of generic families of asymptotically Anti-deSitter soliton-type and black hole 
solutions. 

Besides the numerical results there are also a few analytical ones. Unfortu- 
nately the well developed general theory of dynamical systems, although very 
useful for local questions does not provide the tools required to answer global 
questions in our case. Thus special properties of the equations have to be 
exploited. For negative cosmological constant Winstanley [10] has given an ex- 
istence proof for asymptotically Anti-deSitter black holes with sufficiently large 
radius rh- On the other hand for positive cosmological constant Linden [13] has 
proven the existence of smooth solutions with a regular origin for sufficiently 
small values of A. Clearly it would be desirable to extend Linden's existence 
proof to the numerically found families up to the 3-sphere solutions and to 
black holes, as well as Winstanley's proof to all asymptotically Anti-deSitter 
solutions. 

Experience with the existence proofs for A = [2, 3] suggests that it is quite 
important to classify the possible global behaviour of solutions with a regular 
origin or a black hole horizon. Such a classification is one of the fundamental 
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problems in dynamical systems, and no general method is known to solve it. 

The main result of this paper is a complete classification of all static, spher- 
ically symmetric solutions of the EYM equations with a positive cosmological 
constant. One should also mention here that a partial result in this direction 
has been already obtained by Linden [14] . Since the field equations correspond 
to a rather complicated dynamical system, it is highly remarkable that a com- 
plete global classification of the solutions is possible at all. It should, however, 
be remarked that the numerical results show that the structure of the phase 
space as characterized by the dependence of the solutions on the initial con- 
ditions is rather complicated. This is in particular true inside black hole or 
outside cosmological horizons^ as shown in [15, 16]. 

Our classification proceeds in two steps. We first extend solutions of the ra- 
dial field equations to their maximal interval of existence. In a second step we 
determine the Carter-Penrose diagrams of all 4-dimensional space-times con- 
structible from such radial pieces. 

It may be appropriate to stress the importance of an interplay between 
analytical and numerical methods. On the one hand it seems rather hopeless 
to extract sufficient information on the complex structure of the phase space 
just from studying the field equations analytically. On the other hand precise 
and reliable numerics in particular concerning 'critical' solutions is impossible 
without some analytical understanding of the behaviour of the solutions near 
singular points. 

In Section 2 we set up the field equations. It turns out that through the 
introduction of an auxiliary variable they take the form of a system of coupled 
Riccati equations. In Section 3 we study all singular points and relate them to 
fixed points (f.p.s) of the dynamical system. This requires the introduction of 
different independent and dependent variables in each case, allowing to prove 
the existence and uniqueness of local solutions. In Section 4 we show that any 
solution starting at an arbitrary regular point and integrated in both directions 
either ends at one of the singular points described in Section 3 or runs towards 
r = performing a kind of quasi-periodic motion around a repulsive focal 
point as described in [15, 16]. In Section 5 we show how to glue together 2d 
space-times corresponding to solutions of Eqs. (8) interpolating between two 
singular points and to construct their Carter-Penrose diagrams describing 4d 
inextensible (geodesically complete or singular) space-times. 

Finally, Section 6 contains numerical results illustrating various types of 
solutions. In particular we illustrate how all cases of our classification fit into 
a complete phase space picture. 

The Appendix contains a proposition stating local existence and uniqueness 
of solutions of a class of dynamical systems, allowing to parametrize the solu- 
tions in terms of values at f.p. In addition we formulate several Lemmata used 
in the proof of the classification theorems. 

'^i.e. in the region where the KiUing time is a space-Uke variable 
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2 Preliminaries 



The line element of a spherically symmetric space-time can be written as 

= dsl - r^dn'^ , (1) 

where ds^ is the metric on the 2d orbit space factorizing out the action of the 
rotation group and dJl^ the invariant line element of S'^. The 2d metric can 
always be brought to the diagonal form 

dsl = ai^e^'^dt^ - e^^dB?) , (2) 

where a = ±1 refers to regions, where i is a time-like resp. space-like coordinate. 
Since we consider time-independent solutions it is natural to choose t to be the 
Killing time. This means that we can perform a further dimensional reduction 
to a Id theory. The metric functions z^, A and r then only depend on the radial 
coordinate R. The latter can be chosen arbitrarily exhibiting A as a gauge 
degree of freedom. 

For the SU(2) Yang-Mills field we use the standard minimal spherically 
symmetric (purely 'magnetic') ansatz 

W^Tadx^" = W{R) {Tide + T2 sin 0dip) + Tg cos Odip , (3) 

where Ta denote the generators of SU(2). The reduced EYM action including 
a cosmological constant A can be expressed (in suitable units) as 

S = -J dRe^''+^^ (a + e-^\{r'f + u\r'')')) 

In order to write the Euler-Lagrange equations in first order form we introduce 

N = e-^r', K = re-^v' + N, U = e'^W' , (5) 

and obtain the equations of motion 

re-^N' = {k- N)N - 2U^ , (6a) 
re-^K' = a{l - 2Ar^) + 2U^ - , (6b) 
re'^U' = aWT + {N - k)U , (6c) 

where 

T=- -. (7) 

r 

In addition we have the only remaining diffeomorphism constraint Eq. (9) below. 

We still have to fix the gauge, i.e. choose a radial coordinate. A simple 
choice is to use the geometrical radius r for R (Schwarzschild coordinates), 
yielding e~'^ = N. Putting /x = aN"^, A = e'^/N, and A = Eqs. (6) become 
equal to Eqs. (6) in [3]. This coordinate choice has the disadvantage that the 
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equations become singular at stationary points of r. This problem is avoided 
using the gauge = r introduced in [3], which will also be used here. In 
order to stress the dynamical systems character of Eqs. (5,6) we denote the 
corresponding radial coordinate by r and r derivatives by a dot. We thus get 

r = rN , (8a) 

W = rU, (8b) 

N={k- N)N -2U^ , (8c) 

k = a{l-2Ar^) +2U^ - , (8d) 

U = aWT+{N-K)U , (8e) 

and the constraint (invariant hypersurface) 

2kN -N"^ = + a(l - - Ar^) . (9) 

The metric function v obeys the equation 

v = K-N , (10) 

that can be solved once the other functions are known. Using the constraint 
Eq. (9), the equation for N can also be written as 

iV = |(l-Ar2-r2)-iAf2_f^2_ 

For A = and a = +1 these equations coincide with the Eqs. (49,50) of [3]. 
The quantity T, introduced in Eq. (7) as an abbreviation, obeys 

f = 2UW - NT . (12) 

Viewing T as an additional dependent variable we obtain a system of coupled 
Riccati equations (8) and (12) with two algebraic constraints. This could ex- 
plain, why solutions of this system stay bounded except for simple poles at 
finite r- values (as is true for a single Riccati equation). 
An important quantity is the 'mass function' 

(13a) 



m - 



obeying the equation 



m=(aU^ + —jrN . (13b) 



Furthermore there is a kind of r-dependent 'energy' for the dynamical sys- 
tem Eqs. (8) 

2 

E = ^(2kN - -a{l- Ar^ j^ (14a) 
= ^ - , (14b) 

obeying 
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(14c) 



which wiU prove useful 'Lyapunov Function'. 

Due to the presence of the cosmological constant the vacuum solution ob- 
tained for W'^ = 1 is no longer Minkowski space but deSitter space 

ds^ = fidt^- — -r^dn\ (15) 
A* 

with 

Mds(r) = l-^. (16) 

Because the radial direction changes from space-like to time-like at the cosmo- 
logical horizon Vc = -^3/A we have to subdivide the radial domain into the 
intervals < r < Tc and Tc < r < oo if we use the variables of Eqs. (8). The 
required sign change of a going from one domain to the other can be achieved 
by the substitution (dr, N, k, U) (—idr, iN, in, ill). 

The counterpart of the Schwarzschild solution for A > is the Schwarz- 
schild-deSitter (SdS) solution Eq. (15) with 

, , ^ 2M 

Msds r =1 — , 17 

r 6 

describing a black hole inside a cosmological horizon^ as long as //sas has two 
positive zeros, i.e. M < l/3-\/A. For M = I/Sa/A the function ^g^g has a 
double zero at r = 1/\/A- Geometrically this is the position of an extremal 
(degenerate) cosmological horizon. It corresponds to the fixed point of Eqs. (8) 
{W, U, r, N, k) = (±1, 0, 1/\/A, 0, ±1). As for the extremal Reissner-Nordstr0m 
solution the horizon is at infinite geodesic distance in the space-like direction. 

Likewise there is the Reissner-Nordstr0m-deSitter (RNdS) solution for W = 
with 

2M 1 Ar^ 
MaN(r) = l- — + ^- — , (18) 

carrying an abelian magnetic charge. For A > 1/4 the function fi^^ has only 
one positive zero, while for A < 1/4 it has three positive zeros as long as 
M_(A) < M < M+(A) for certain functions M±(A). For M = M_(A) the 
two smaller zeros merge to a minimum of fi^jq at r?. = (1 - V1-4A)/2A, 
while for M = (A) the two larger ones merge to a maximum of ^Urn at r\ = 
(1-FV1 - 4A)/2A. For A = 1/4 and M = 2a/2/3 ah three zeros merge to a triple 
zero at r = \/2. As long as ji^m has three positive zeros the solution describes 
a charged black hole located inside a cosmological horizon. The horizon at 
the smallest value of r, a cosmological horizon according to our terminology, 
corresponds to the Cauchy horizon of the Reissner-Nordstr0m black hole. The 
two possible double zeros at r = r± are the positions of extremal horizons of 
the black hole resp. cosmological type. Only the extremal black hole horizon 
has a limit for A — ^ 0. 



^We call a horizon where r grows, while a in Eq. (2) becomes positive a black hole horizon 
and a horizon, where r grows, while a becomes negative, a cosmological horizon. 
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Apart from these general solutions of deSitter type there are some degener- 
ate ones with r and W constant (and henee N = U = 0) and a space-time of 
the form M2 x S"^. In fact these solutions correspond to double zeros of //g^s or 
/XaN, resp. to the SdS and RN fixed points (f.p.s) discussed in Section 3.8. Prom 
Eqs. (8,9) one finds that either = 1 (SdS) or = (RN) and k = - 
with 

= a(l - 2Ar2) , and 2Ar2 = P ^ ^ for ^ 1 , ^ 

^ ^ ^ ' \l± VI -4A for = . ^ ^ 



For > there are two solutions 

k{t) = Kg tanh(Ksr) , 
k{t) = Kg coth(Ksr) , 

with the hne element ds^ = a{e^^dt^ 
[9, 17]. For k1 = —c? < one obtains 

k{t) = a cot (ar) , 



e'^" = cosh2(ACsr) , (20a) 
e'^" = K-'^smh'^{KsT) , (20b) 

r'^dr^) known as the Nariai solutions 
e^^^ = a-2sin2(ar) , (20c) 



and for = in the limit A ^ 1/4, ^ 

ac(t) = i , e^^^ = (20d) 

Geometrically the first three are SO {2, 1)/S'0(1, 1), while the last one is flat 
space (in unusual coordinates). 

3 Classification of Singular Points 

Our aim is to give a complete classification of all solutions starting at an arbi- 
trary regular point. We may integrate Eqs. (8) in both directions until cither 
one of the dependent variables diverges for finite r or r ^ ±00 where the solu- 
tion may converge (to a finite or infinite value) or not. As the experience with 
a single Riccati equation shows we have to expect singularities for finite values 
of r. In fact, a typical example is given by non-degenerate horizons, where the 
function k has a simple pole. Natural candidates for the convergent case for 
T — > ±00 are solutions running into one of the fixed points of Eqs. (8). The 
latter are easily found due to the simple structure of these equations. With 
exception of one case discussed in Section 3.5, the singularities occurring at 
finite values of r turn out to be simple poles. These singularities can all be 
transformed into f.p.s by suitable changes of the independent and dependent 
variables. Furthermore all the f.p.s except the RN f.p. for A = 1/4 are hyper- 
bolic ones. Using a technique known as 'blow-up' [20, 21] we are able to reduce 
also this exceptional case to the hyperbolic one. 

Reducing the study of singular points to that of hyperbolic f.p.s has the 
advantage that the theory of dynamical systems provides powerful theorems 
about the local behaviour of solutions in their neighbourhood [18, 19, 20]. The 
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theorem of Hartman and Grobman about the existence of a stable (and unsta- 
ble) manifold guarantees local existence and uniqueness of solutions. In order 
to parametrize the solutions by values determined at the f.p. itself, we rewrite 
the equations in the (non-autonomous and non-hyperbolic) form required by 
Prop. 12 in the appendix, whenever possible. 

Most of the singularities are actually of geometrical origin. Since our in- 
dependent variable r parametrizes the space of orbits of the isometry group 
G = i? X 50(3) of time-translations and rotations acting on the 4d space-time 
M4, we expect singularities at the fixed points of the group action. For rotations 
this is r = 0, the center of symmetry. It turns out that the point r = can be 
reached for finite r, if |A^| 00 (curvature singularity) or for r — > ±00 with 
iV — ^ ^1 (regular origin). Regular (non-degenerate) horizons are represented 
in the quotient space by their 'bifurcation surface' [22] consisting of f.p.s of the 
time-translations. The required vanishing of e'^ is possible only if k diverges at 
some finite r. In contrast, the degenerate horizons we find require r — ±00 
and finite k. A deSitter type singularity reached for finite r occurs for r — >■ 00 
in the region a = —1 due to the presence of the cosmological constant. 

In the subsequent paragraphs we shall give a characterization of the so- 
lutions near the singular points discussed above, providing the basis for our 
classification theorems. 



3.1 Regular Origin 

As a first case there are solutions for cr = +1 with a regular origin r ^ 0, 
W"^ 1, and A'' ±1 for r ibcx). These are f.p.s of Eqs. (8,9) with two 
dimensional stable manifolds. In order to characterize the corresponding one 
parameter families of solutions by a value h determined at the f.p., one rewrites 
Eqs. (8,9) in the form required by Prop. 12 as in [3]. Choosing W ^ +1 without 
restriction one thus obtains unique solutions which are analytic in h and r with 

W{r) = l-hr'^ + O(r^) , (21a) 
Ar2(r) = l- (^452 + I) r2 + 0(/) . (21b) 

3.2 Reissner-Nordstr0m Singuleirity at r = 

The generic behaviour at r = is singular in complete analogy to the case 
A = [3], since the terms proportional to A may be neglected for r — 0. 
The behaviour near r = depends on the sign of a. For a = +1 one finds a 
Reissner-Nordstr0m type singularity of the metric 

W(r) = W„- + ly-a'-" + 0{r*) , (22a) 

iVV) = 2^^°-^ + 0(l). (22b) 

with three parameters Mq, Wq, and W3. As was shown in [16] the singularity 
may be described as a f.p. introducing a suitable set of dependent variables. 
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Here we shall give an alternative formulation using Schwarzschild coordinates 
and the dependent variables already introduced in Prop. 13 of [3]. 

r^N = r (kN - 2rU'^) /N , (23a) 



dr 

r-^K = r(l-2Ar^ + 2U^ - K^^/N , (23b) 

r-^W = r^U/N , (23c) 
dr 

r^U = r{\- kU^/N , (23d) 

r^X = r{^W'^ -2U'^ -l^U/N , (23e) 

with N = rN and A = WT + NU . These equations are precisely of the form 
required for the application of Prop. 12 implying the regularity of local solutions 
as a function of r obeying the constraints 

iV(0) = ±(H"2^0) - 1) , iV(0)C/(0) = W^(0)(W^2(0) - 1) , (24) 

following from the finiteness of A(0) and the constraint Eq. (9). Geometri- 
cally this behaviour leads to a Rcissncr-Nordstr0m type curvature singularity 
at which the space-time is geodesically incomplete. 

3.3 Pseudo-Reissner-Nordstr0m Singularity at r = 

There is an analogous singularity for cr = — 1 called pseudo-RN singularity in 
[16], with a rather different behaviour of the metric and of the YM potential 
W compared to the case a = +1. 

W{r) = Wo ± r + 0(r2) , (25a) 

„V)=(lIl_i£^l!Wil^^o(i), (25b) 

The description as a f.p. was already given in [16], but here we give an alterna- 
tive description, providing a unified treatment of the two different signs of W 
at r = 0. We introduce the dependent variables U = U/N, T = T/N, 

[72-1 2WUf'^ T2-1 

^ = —;— + ' ^^^^ 



Using Schwarzschild coordinates one finds 







r-^W = rfw , 
ar 


fw 




fu 


r^U = rfu , 
ar 




h 


r-^f = rfr , 
dr 



U , (27a) 
+ 2zU , (27b) 



W'^ - 1 

2WU{T^-l)-\ 



(27c) 
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r-^x = -x-2y + rf^c , fx = -2xy + 2 [2(x + y + 2z)WU + 
ar L 



(1 - Ar2)c72 + W^{2U^ _ r2 - 4U^f^) 
r—y = 2x-2y + rfy, fy = 2{x + z)y 



f2 



+ 



rl - Ar 



2 



AWUy 



(27d) 



(27e) 



with 2; = [r(l— Ar^)T'^/(VK^ — 1)^— y]/2. In their linearized form these equations 
have zero eigenvalues for W, U, and T, while x and y contribute with the 
complex conjugate eigenvalues (—3 it ■\/l5i)/2 representing divergent modes. 
Introducing the linear combinations (eigenmodes) x± = {1 i: ■\/TEi)x — 4y we 
obtain 

r^x±= ~^^^^ x± + rU, U = {l±Vl5i)U-4fy, (27f) 

satisfying the requirements of Prop. 12. We thus obtain local solutions of 
Eqs. (27) with x(0) = y(0) = depending analytically on W{0) ^ ±1, C7(0), 
T(0), and r. Imposing Eqs. (26) as constraints we then obtain the conditions 
[72(0) = f 2(0) = 1 and from the constraint Eq. (9) we find k/N = 1+U'^-rz 
2 for r ^ 0. 

3.4 Schwarzschild Singuleirity at r = 

The form of the RN type singularities requires that W^{0) 7^ 1, while for 
W'^{0) = 1 we obtain a Schwarzschild type behaviour {N'^ ~ 1/r) with the 
sign of the mass determined by a. This case requires a different linearization, 
which we shall give subsequently. Assuming without restriction W{0) = +1, 
we introduce the dependent variables W = {W — l)/r^, U = U/rN, k = k/N, 
N = avN"^, and y = {\ — 2k/N + 2r'^tj'^)/r, use the constraint Eq. (9) to express 
y = iV-i(l - r2(2 + r'^WfW^ - Ar^), and find 

r^N = rU, U = {y-2rU^)N, (28a) 

r-^W = U -2W , (28b) 
ar 

ar N 

In order to bring Eq. (28b) into the form required by Prop. 12, we further 
introduce W = W - U /2 with 

r-^W = -2W -'^^ , (28d) 

and obtain solutions with {W,tJ,N) tending to finite values {—b,—2b,—2M) 
with aM < 0. The corresponding asymptotic behaviour at r = is 

W = l-br^ + O(r^) , (29a) 
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aN'^ = + l + 0(r), (29b) 

'^ = y (1 + 0(0)- (29c) 

The RN and S type singularities at r = are different singular points witfi 
different asymptotic behaviour. The S case (for a = +1) can, however, be 
seen as the limit Wq ^ ±1 of the more general RN case. In Fig. 1 we present 
numerical solutions for A = 0.8 starting at a regular origin with h ~ 0.237 as 
they approach the RN type singularity with Wq k. —1. N and k grow like l/^/r 
and i7 ~ as for the S type singularity in an interval that increases as \Wq — 1| 
decreases until eventually N starts to grow like 1/r whereas k has a finite limit 
and U tends to =bVFo- Thus in the limit \Wq\ — > 1 the functions W{r), U{r), 
N{r), and K{r) converge pointwise to those for the S type singularity. 



3.5 Unbounded Oscillations at r = 

The nature of the f.p. is rather different in the RN and pseudo-RN case. While 
for (J = +1 the f.p. is a simple attractor it is a partially repulsive (in the U, T 
subsystem) focal point for a = —1. 

The RN singularity at the origin for cr = +1 is therefore generic (i.e. persis- 
tent for sufficiently small changes in the initial data and A), as is the deSitter 
asymptotics for a = —1 described below. All other singular points discussed 
above and below are not generic because they require one or more divergent 
modes to be suppressed. 

The generic behaviour for a = —1 and r decreasing is somewhat more 
complicated. As described in [15, 16] this leads to an oscillatory behaviour of 
the solutions as they approach r = with k — > — oo and unbounded oscillating 
U and N. We introduce new dependent variables n = N/k, u = U/k, and 
t = T/k which are bounded due to Eq. (9) since k diverges. In addition we 
introduce a new independent variable p with dp = \K\dT such that p ^ oo as 
K —>■ —DO (see the proof of Prop. 10 for details) and denote the derivative with 
respect to p by a prime. Using z = 1/|k| we find 

r' = —rn , (30a) 

u' = iylv? - n + (2Ar2 - \)z^) u - Wzt , (30b) 

t' = {2u^ + n - 1 + (2Ar2 - l)z^) t + 2Wzu , (30c) 

n' = (2t(2 + n - 2 + {2Ar^ - l)z'^) n + 2u^ , (30d) 

z' = (2u2 - 1 + (2Ar^ - l)z^) z , (30e) 



and the constraint 



2n = n^ + 2u^ + t^ + z'^ (Ar^ - 1) , (30f ) 

with \W\^/z = \J z — rt bounded. Neglecting terms that vanish as z — > we 
obtain the asymptotic equations 

u' = (2«2 -n)u, (31a) 



11 



t' = {2u'^ +n-l)t, (31b) 

n' = (2u^ +n-2)n + 2u'^ , (31c) 

with the constraint 

2n = + 2u^ + , (31d) 

or, in terms of new dependent variables ^ = n = N/k, rj = 2n^/n = 21/"^ /Nk, 
and C = t'^/n = T^/Nk 

| = e-2 + (e + l)^ = 4(2-0-(C + l)C, (32a) 

!L = 2-n + {n-3)c = -{2-nf + {3-n)c, (32b) 

V 

| = C-2 + (4-^-C)e = 2-C-(^ + C)^, (32c) 
with the constraint 

2 = C + V + C- (32d) 

Fig. 2 shows four trajectories as they emerge from the focal point {^,r),Q = 
(1/2, 1, 1/2) corresponding to the pseudo-RN singularity and rapidly approach 
the bounding triangle consisting of the three separatrices ^ = 0, ?? = 0, and 
C = connecting the three hyperbolic fixed points A = (2, 0, 0) corresponding 
to the S singularity, B = (0, 0, 2) corresponding to a regular horizon described 
below, and C = (0,2,0). 

The behaviour referred to as 'chaotic' in [16] is due to this rapid approach to 
the boundary. The terms proportional to W in the full Eqs. (30) may become 
relevant when rj <^ z or <^ z. A solution may then either end with an S 
singularity or regular horizon or nearly miss them with a resulting substantial 
change in t or u. There is no such instability for ^ <C 2. 

3.6 Regular Horizons 

In the case of a regular, i.e. non-degenerate horizon we may use s = 1/k as the 
independent variable tending to zero. We thus replace Eqs. (8) by 

s^T = -'^ = s{l + s'^F{r,U)) , (33a) 

CIS K 

s^r = s{l + s^F{r,U))rN , (33b) 

s^W = s(l + s^F(r,U))rU , (33c) 
as 

s^U = -U + s{l + s^Fir, U)) iaWT + NU) - s^FU , (33d) 
ds 



with 



^) = r^:^ ' / = ^ (1 - 2Ar2) + 2C/2 , (34) 
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while the function A'^ can be computed from the constraint Eq. (9). These 
equations have precisely the form required in Prop. 12 with C/ — > and r, r, 
W tending to finite values r/j, r/j, W^. 

In terms of the coordinate r the behaviour near the horizon is (performing 
a shift in r we arrange for = 0) 

r{r) = rh(^l + N, y ^ + 0(t") , (35a) 

N{t) = iVi(r + 0(t^)) - Wfr^ + O(r^) , (35b) 

Wir) = Wh + rhWi y + O(t^) , (35c) 

where 

iVi = ^(l-Ar2-r2)| , Wi = ^WT . (36) 

Z \h Z h 

According to our definitions we have a black hole horizon if crA'^i > or a 
cosmological horizon if aNi < 0. If Ni vanishes the expansion of N starts with 
the cubic term with a corresponding change in the expansion of r. Geometrically 
this describes the situation that the horizon coincides with an equator, i.e. a 
maximum of r. The condition o'A'^i > determines the domains C-). = {1 — 
Ar^ — > 0} in the (r, W) plane. In terms of and the boundary curves 
are (parts of) ellipses. In Fig. 3 these domains are shown for various values of 
A, for A > 1/4 they become disconnected. 

Solutions with a regular origin can be seen as a limit r/i ^ of those with 
a regular horizon. In Fig. 4 we present solutions for A ~ 0.3642 starting at 
a regular origin with b ~ 0.4296 (i.e., near the point i?2 of Fig. 14) as they 
approach a regular horizon of decreasing radius. ~ k ~ — 1 and U/r ~ 26 as 
for the regular 3-sphere solution with 2 zeros of W in an interval that increases 
as r/i — until eventually k starts to diverge whereas N and U tend to zero. 
Thus the black hole solutions converge to the one with a regular origin. 

3.7 deSitter Asymptotics 

For (7 = — 1 the radial variable r can tend to oo leading to deSitter asymptotics 

W = Wo + ^+0{l/r^) , (37a) 
N^ = ^r^^l + ^ + 0{l/r^) . (37b) 

We use s = l/r as independent variable and introduce the dependent variables 
U = UN, N = N/r, and the 'mass function' Eq. (13). 

'is^='i^^^'-^^+'^''w) ' ^^^^^ 

with iV^ = A/3 — + 2s^m. Thus all the requirements of Prop. 12 are fulfilled. 
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3.8 The Schwcirzschild-deSitter and Reissner-Nordstr0m 
Fixed Points 



Besides the singular points reached at finite r there are the f.p.s of Eqs. (8) with 
{r,W,U, N, n) (rs,Ws,Us, Ng, Kg) as r ^ oo. One such f.p. is the regular 
origin = already discussed. For finite / this requires Ug = Ng = and 
either Wg = 1 or Wg = 0, which we denote as Schwarzschild-deSitter (SdS) resp. 
Reissner-Nordstr0m (RN) f.p. The finite values of Kg and rg are determined by 
Eqs. (8d) and (9). In order to prove the local existence of solutions running 
into these f.p.s we will apply the usual linearization technique, using variables 
f = r/rg — 1, W = {W — Wg)/r, and R = k — Kg. 

The SdS f.p. requires cr = — 1 in order to get > and we find rg = l/s/A 
and Kg = ±1. Prom Eqs. (8) we obtain 

N = KgN + fi, (39a) 

R = -2Kg{K + N) + f2, (39b) 

W = U + f3, (39c) 

U = -2W - KgU + /4 , (39d) 

and the constraint Eq. (9) yields f = N/Kg + f^. The functions z = 1, . . . , 5 

are analytic and at least quadratic in the dependent variables W, U, N, and 
R. The eigenvalues of the linear system are Kg = ±1, —2ks = =f2 for the N, 
K part and {—Kg ± \/7«)/2 in the W, U sector. For Kg = 1 there are three 
convergent modes for r — ^ oo, while for Kg = —1 there is only one convergent 
mode. For Kg = 1 the function W performs damped oscillations of period 
47r/\/7 around the limiting value Wg = ±1. For Kg = —1 both divergent 
modes of the W, U subsystem have to be suppressed. However, the r.h.s. of 
Eqs. (39c,d) vanish fov W = U = and thus the stable manifold of the f.p. 
contains only the SdS solution Eq. (17) with W = ±1 and a double zero of 
Hsds- Similarly for Kg = +1 the divergent mode of N has to be suppressed and 
thus N{t) « /~ 2U'^e^-^'dT' > 0. Except for the Nariai solutions (20a,b) with 
W = f = 0, the f.p. can only be reached from r < rg with N > 0. 

The RN f.p. with Wg = leads to = = (1 ± V1-4A)/2A and 
= (j(l — 2Ar^). It requires A < 1/4 and the combinations a = +1, r_ resp. 
(7 = — 1, r+ to obtain k"^ > 0. For the linearization we assume A < 1/4, i.e. 
Kg ^ and obtain 

N = KgN + h, (40a) 
k = -4(7 Ar^— - 2KgR + /2 , (40b) 

Kg 

W = U + f3, (40c) 
if = -aW - KgU + fi , (40d) 

where f = N/Kg + is again obtained from the constraint. The eigenvalues of 
the linear system are Kg, —2Kg for AT, k and {—Kg it \/ k1 — 4cr)/2 for the W, U 
system. Since k^ < 1 the latter eigenvalues are complex conjugate for a = +1 
and real for cr = —1. For cr = +1 there are three convergent modes for Kg > Q 
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and only one convergent mode for Kg < 0, whereas for a = —1 there are always 
two convergent modes. The need to suppress the divergent modes for a = +1 
and Kg < again allows only the local (and thus global) solution W = 0, i.e. the 
stable manifold contains only the RNdS solution. The period of the damped 
oscillations of W around zero in the case a = +1, Kg > is ^tk j ^ A: — and, 
as in the SdS case, the f.p. can only be reached from r < rg with A'^ > 0, except 
for the Nariai solutions with W = f = 0. For a = —1, Kg > there are again 
the Nariai solutions and solutions with W ^ 0, N > 0, and r < rg. For a = —1, 
Kg <0 there is the RNdS solution and there are solutions with W ^ 0. 

3.9 The Reissner-Nordstr0m Fixed Point for A = 1/4 

For A = 1/4 some of the eigenvalues of the linearized Eqs. (40) vanish and 
this degenerate case requires special consideration. First we introduce new 
dependent variables 



yy, = W = W/r, yu^U = U + KW/2, 
y^ = K, y, = f=l-rV2, yM = N = N + WU ^ > 



and in addition 

yA = A = aW^ + . (42) 
Next we rewrite Eqs. (8) and the constraint Eq. (9) as 

yi = Fiiy) + My) , i = W,U,K,r,N,A , (43a) 
= Fe(y) + /e(y) , (43b) 

with the 'leading' contributions 







Fu) 




Fn 


= af — 


Fr 


= -2N , 


F^ 




Fa 


= -nA , 


Fc 


= A-kN 



1 \(w 

-k/2) \ U 



(44a) 

(44b) 
(44c) 

(44d) 
(44e) 

4 ' (44f) 
dominating the behaviour near the f.p., and with the 'non-leading' corrections 

(-N 0\ (W\ J ar 



fuJ^\fuNl\Uj' fu = ^-^+<yr'W' + U'-KN , (45a) 

/k = 2?72, (45b) 
fr = 2{WU + fN) , (45c) 



ar 



3 iV2 ar^W^ 



fA = 2N{U'' -aW'') + 2fuWU , (45e) 
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Prom the linearized form of tlicsc equations we get vanishing eigenvalues from k, 
f, N, and A, whereas W and U contribute the eigenvalues zti for a = +1 or ±1 
for £7 = —1. We therefore expect that (at least some of) the functions decrease 
with powers of r instead of exponentially as for non-vanishing eigenvalues. 

In order to analyze the behaviour of solutions tending to the f.p. as r — ^ oo 
we perform a 'quasi-homogeneous blow-up', following [21] with some minor 
variations. In terms of rescaled variables 

yi = A~"'yi , a„ = 1 , = Ojj = ar = 2 , cnjv = 3 , = 4 , (46) 

we find 

F,{y) = X''i+'Fi{y) , /,(y) = 0(A"'+2) , for i = K,r,N,A, (47a) 

F,(y) = A^Fe(y), f,{y) = 0{X'). (47b) 

Furthermore the leading contributions Fi, i = K,,r, N, A and F^. depend only on 
K, f, N, and A. Hence they define a quasi-homogeneous dynamical system for 
{k, f, N, A) of type a = (1, 2, 3, 4) and degree 2 in the sense of [21], whereas the 
non-leading corrections yield at least one additional power of A upon rescaling. 
The non-leading correction terms, fi depend on W and U individually, not only 
on the combination A. One can easily find bounds for \W\ and \U\: for cr = -|-1 

they are bounded by whereas for cr = — 1 with eigenvalues ±1 they are 
bounded by exp(— 5t) for any S < 1 due to the special form of Eqs. (44a, 45a). 

Next we introduce jii = Mjai with some M such that /Xj > 2, new dependent 
variables A and y 

= E 7i \yiT' , Vi = y^-'^'Vi , i = K,r,N,A, (48) 

i 

with some constants 7i > such that J2i 7i > Oj ^ i^^w independent variable 
t with dt = Xdr. For 7i > this is a 'polar' blow-up whereas e.g., 7 = (1, 0, 0, 0) 
yields a 'directional' blow-up in the K-direction. Neglecting for a moment the 
non-leading correction terms we obtain 

1 dX sr^ , ,1,- FAy) , , 

T:77 = E7^ y^r^' 49a 
A dt ^ aiyi 

_| = F,(,) (49b) 

i.e., an y-dependent 'radial' flow with the f.p. A = and a A- independent 
'tangential' flow with f.p.s determined by the conditions 

amFjiy) = ajyjFi{y) , i,j = k, r,N,A . (50) 

There are three pairs of tangential f.p.s Pi±, i = 1, 2, 3 with k^O and 

Pi± : r = N = A = , (51a) 

„ ~ 3(7 _9 ~ 3(7 _Q r 21a ^ A , 

P2±: r = —i^, N = —k\ A = —k\ (51b) 
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Inserting these values into Eq. (49a) yields dX/dt = —kX at Pi±, — kA/4 at 
P2±, and +kA/3 at P3±. Thus Pi+, P2+, and P^^ arc attractors for the radial 
flow, whereas Pi-, P2-, and P3_|_ are repulsive. On the 2-dimcnsional surface 
y subject to the constraint Fc{y) = the f.p.s Pi- and P3+ are attractors 
for the tangential flow, P^. and Pa- are repulsive, and P2± are saddle points. 
Furthermore there are trajectories along the curve N"^ = ar^ jVl connecting all 
six f.p.s as shown in Fig. 5: separatrices from Pi_|_ and Pa- to P2+, separatrices 
from P2- to Pi- and P3+, as well as trajectories from Pi_|_ to P3+ and from Pa- 
to Pi-. This leaves no room for limit or separatrix cycles, therefore according 
to the Poincare-Bendixson theory the f.p.s are the only limit points of the 
tangential flow. 

With this blow-up we have transformed the degenerate f.p. y = into a set 
of hyperbolic f.p.s (A, y) = (0, Pj±) with one dimensional stable manifolds for 
(0, Pi+), (0,P2-), (OjPa-) and two dimensional stable manifolds for (0, Pi_), 
(0, P2+), (0, Pa+). The non-leading correction terms, /j in Eqs. (43) yield 0(A) 
corrections to the r.h.s. of Eqs. (49) and thus preserve the position and character 
of the six f.p.s. The structural stability of the system (43) guarantees that the 
flow pattern is not changed by the non-leading corrections. Moreover the stable 
manifolds for (0, Pi-), (0, P2-), and (0, Pa+) are purely tangential at A = and 
therefore excluded for trajectories tending to the f.p. y = as r — > 00. 

We have chosen W,U, and N such that l/X — Cikr has a finite limit on each 
of the stable manifolds with ci = 1 for Pi+, C2 = 1/4 for P2+, and ca = —1/3 
for Pa_ . Consequently we choose yet another independent variable s = A = 1/r 
as well as dependent variables y as in Eq. (46) to obtain 

d_l'W\^fK/2-2-l/s\fW\_ fgwO\fW 
^ds\Uj~\ a/s R/2-2)\u) ' \gu 9u ) \U 

s^K = k^ -k-ar- sg^ , g^ = l^^ (52b) 

ds s"^ 



f-^r = 2{N -r)-sgr , ^ ^ 



s^^N = ^-SN-sg^ , gN = ^, (52d) 

s^A ={k- 4)A -sg^, = 7^ ' ^^^''^ 

as well as the constraint 



= A-kN- — + sgc, gc=^, (52f) 



with g^s that are analytic in s and the rescaled functions. The requirement 
that the r.h.s. of Eqs. (52b-e) vanish at s = yields the conditions Eqs. (51) 
for the f.p.s with k{0) = 1 for Pi+, k{0) = 4 for P2+, and k{0) = -3 for Pa_. 
These three different solutions are to be expected as the limit A 1/4 of the 
Nariai solution (20b), and the RN f.p.s with the combinations r+, k > and 
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r_, K < 0. In terms of shifted variables r = f — f(0) etc., that vanish at s = 
we thus obtain the hnearized equations 
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(53) 



with eigenvalues +1 (corresponding to a shift in r), /t(0) — 4 (due to the con- 
straint), as well as —2, —3 for Pi+, +7, —6 for P2+1 and —6, —7 for Ps-. 

We still need an additional equation for the YM functions W and/or U . For 
(7 = +1 the solutions with k{^) < 4 yield W = whereas k(0) = 4 requires 
W ^ 0. Motivated by the eigenvalues ±i, we introduce in the latter case 
W = VAsm{T + 6) and U = \/Acos{t + 9). The equation 



Qw) cos(r + 9) + Qu sin(r + 6*)) sin(T + 9) 



(54) 



implies that the phase 9 has a finite limit as s 0. 

For a = —1 with the eigenvalues ±1 we introduce W± 
Eq. (52f) by 



W and replace 



= ^ + ^) W± - s{g±W± + g± 



1 K 

2 



(55) 



where g± and ^± are linear combinations of g-w, Qu-, and gu. The divergent mode 

W+ has to be suppressed and consequently Wj^ and decrease exponentially, 
implying that A = W^W- vanishes at s = and thus excluding the solution 
with k{0) = 4. As for A < 1/4, there are solutions with k > and k < with 
both W = OandW^(} 



4 Classification of Solutions 

Starting from any regular point r = tq of the Eqs. (8) we may integrate these 
equations to both sides until we meet a singular point. In case / r dr is finite the 
corresponding space-time is geodesically incomplete, otherwise it is complete 
in the radial direction. This does not necessarily mean completeness in all 
directions as the example of the extremal RN solution demonstrates. A general 
recipe how to glue together incomplete pieces of space-time to a geodesically 
complete manifold is described in [23]. 

In the preceding section we have already discussed the possible singular 
points we may encounter integrating Eqs. (8) and described the local behaviour 
in their neighbourhood. The missing step is a classification of the singular 
points reachable from some regular point integrating the equations to both 
sides. Instead of integrating towards decreasing r we may as well reverse the 
sign of N, K, and U and always integrate towards increasing r. 

The global classification is different for the two possibilities a = ±1, i.e. for 
regions with space-like or time-like radial direction. Here we start with some 
general properties equally valid for both cases. 
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Proposition 1: Consider a solution of Eqs. (8,9). 

i) If N{to) < for some tq then N{t) < for all r > tq. 
a) r has no minima and at most one maximum. 

Hi) r has a limit ri as t ^ Ti with In ri finite if and only if J^^ N dr is finite. 

iv) If N{t) and k{t) are finite at some r then W{t) and U{t) are finite. 

v) If N < 0, K < 0, and k + N is unbounded, then k + N ^ — oo at some 
finite ri . 

Let us define = {k — N)/{k+N) and 7 = {t — ti){k+ N). IfC, is constant 
then 7 = 4/(1 + 3^^) + 0((r — ti)^). If C, is restricted to an interval then 
7 is some (t dependent) average of this expression over that interval, in 
the most general case 1 < 7 < 4 up to terms 0((r — ti)^). 

Proof: 

i) This is an immediate consequence of Eq. (8c). 

ii) This is an immediate consequence of the previous point and Eq. (8a). 

iii) Property ii) imphes that r is monotonic for sufficiently large r and there- 
fore has a limit; the finiteness of In r is a consequence of Eq. (8a). 

iv) Due to the previous point Inr has a finite limit and U"^ is integrable due 
to Eq. (8c). The Schwarz inequality then implies that \U\ and thus r\U\ 
are integrable, i.e. W is bounded. The boundedness of U finally follows 
from Lemma 13 in the appendix applied to Eq. (Be). 

v) Prom Eqs. (8c,d) we obtain 

{k + Nj = a{l - 2Ar^) -^{K + Nf-^{K- Nf 

<a{l-2Ar^)-^{K + Nf , (56) 

with 1 — 2Ar^ bounded since N < 0. Therefore Lemma 14 implies that 
K+N diverges for some finite ri and we assume ri = without restriction. 
Introducing r/ = 1/{t{k + N)) we find 

ri{r) = J\l + 3C')dT' + ar (1 - 2Ar^WdT' . (57) 

Since < 1 + O(r^) due to Lemma 14, the second integral on the r.h.s. 
is O(r^), whereas the first integral yields I/7 as the average value of 
(1 + 3e)/4. U 

According to Prop. 1 we can choose a suitable fo > tq such that has a 
definite sign for r > tq. We have to consider the cases that either N and k are 
bounded for all r > tq or that A^ and/or n are unbounded. 

For N and k bounded we have the result: 



19 



Proposition 2: Suppose a solution of Eqs. (8,9) with N, k, W, U, and Inr 
bounded for r > tq, then the solution tends to one of the f.p.s of Eqs. (8) as 
T oo: 

i) for a = -1, the SdS fp. described by Eqs. (39); 

ii) for a = ±1, the RN fp. described by Eqs. (40). 

Proof: Since Inr is bounded, |A^| is integrable and tlierefore is integrable 
due to Eq. (8c). Integrating 

(U^y = 3(aWT +{N - k)U^ , (58) 

derived from Eq. (8e), we see that U has a hmit, which must be zero since 
otherwise W would diverge. Lemma 14 applied to Eq. (8d) implies that k has 
a limit. From Eqs. (14) then follows that the 'energy' E is bounded and \E\ 
is integrable, thus E has a limit. This implies that N has a limit which must 
be zero since otherwise Inr would diverge, and that W has a finite limit which 
must be one of the f.p.s of Eq. (8e), — ±1 for the SdS f.p. or W ^ for the 
RN f.p. ■ 

For N bounded and k unbounded we have the result: 
Proposition 3: Suppose a solution of Eqs. (8,9) with N bounded and k un- 
bounded, then the solution has a regular horizon at some finite ti as described 
in Eqs. (35,36). 

Proof: First we show that ri must be finite. For cr = +1 we use Eq. (56), for 
(T = — 1 we use 

{K + Ny=l-AU^-2T^-K^ + 5KN-3N'^ < 1 + 7N{k + N) - {k + Nf . (59) 

In both cases Lemma 15 implies that k + N and thus k diverges to — oo for 
some finite ri, thus Inr has a finite limit. We assume again ri = without 
restriction. 

Next, for a = +1 Eq. (11) implies that U'^ is integrable, since otherwise 
N would be unbounded from below. Therefore \U\ is integrable due to the 
Schwarz inequality and thus W has a finite limit. The constraint Eq. (9) then 
implies that U^/k is bounded. For a = —1 the constraint implies that U^/k is 
bounded. 

For a = ±1 we use s = 1/k as in Eqs. (33,34) with dr = (l + s'^F) ds and 
find that sf and hence sF is bounded. Therefore \U\dT < oo and thus W 
has a finite limit. We now turn to Eq. (33d) for U and obtain 



^ ^..'^g) + -'> , (60) 

as \ r / ^ ' 



with G{s) = — Jg s''^F ds' bounded, yielding 



(1 + ^"F) Us' . (61) 



r s Jg s r 
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The constant c must vanish, otherwise sU"^ would not be bounded. This imphes 
[7 ^ for ,s ^ 0. From the constraint we finally obtain N(0) = 0, thus the 
solution has a horizon. ■ 

4.1 Solutions with cr = +1 

To begin with we state some properties of the solutions required for their global 
classification. Most of these are trivial generalizations of results already proved 
for A = in [3], but we repeat their proofs for completeness. 

Proposition 4: Consider a solution of Eqs. (8,9) with a = +1. 

i) The function W can have neither maxima if W > 1 or Q > W > —1 nor 
minima if W < —1 or < W < I. 

a) N is bounded from above and if N{to) < 1 for some tq then N{t) < 1 for 
all T > Tq. 

Hi) r is bounded from above. 

iv) K + N is bounded from above and if {k + N){tq) < 2 for some tq then 
{k + N){t) < 2 for all r > tq. 

v) If W'^{to) > 1 and WU(to) > for some tq and k is bounded, then 
N{ti) < —1 for some Rnite t\ > tq. 

vi) If N{to) < — 1 for some tq and k is bounded, then N{t) diverges to — oo 
for some finite t\ > tq. 

Proof: 

i) This is an immediate consequence of Eqs. (8b ,e). 

ii) This follows from Eq. (11). 

iii) In the following argument we can always assume Ar^ > 3 and iV > since 
otherwise there is nothing to prove. Thus < — 1 due to Eq. (11) and 
therefore A^ = for some finite ri where r has a finite maximum. 

iv) Eq. (56) implies (k + ATJ < for k + iV > 2. 

v) In the following argument we can always assume A^ > — 1 since otherwise 
there is nothing to prove. Let ^ > 1 be a constant such that k < 2A. 
First we note that W'^ monotonically increases and WU > for r > tq. 
From W = W{W'^ ~ 1) + (2A?^ — k)W we obtain a positive lower bound 
for 1 1^1 and thus there exists a tq > tq such that |W^|(to) > A. Using the 
Schwarz inequality we find 

\TU\-=2\W\{2T^ + U'^)- k\TU\ > {VS - 2)\TU\ . (62) 

Therefore \TU\ and, using again the Schwarz inequality, 2T^ + grow 
exponentially for r > tq. Eq. (11) then implies that A" < —1 for some 
finite Ti. 
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vi) This is a consequence of Eq. (11) and Lemma 14. ■ 

In addition there are new features compared to the case A = 0, which are 
mainly due to the possibihty that the variable k may become negative. 

Proposition 5: Suppose a solution of Eqs. (8,9) for a = +1 with N and k 
bounded for all r > tq, then the solution runs into 

i) the RN f.p. described by Eqs. (40), if r ^ ri > (with Kg < only- 
possible ifW = ); 

a) a regular origin as described by Eqs. (21 ), if r ^ 0. 

Proof: First we note that the limit ri of r cannot exceed 1 / a/A since otherwise 
N would diverge due to Eq. (11). Moreover the 'mass function' Eqs. (13) is 
monotonic and bounded and therefore has a limit. Furthermore W is bounded 
due to Prop. 4. 

i) If ri > then U is bounded due to the constraint Eq. (9) and Prop. 2 
implies that the solution runs into the RN f.p., as discussed in Section 3. 
For reasons given there the case /c^ < is only possible ifW = 0. 

ii) For ri = we first of all choose a tq > tq such that A'' < and Ar^ < 1/2 
for all r > To and restrict the following arguments to this subinterval. 

Eq. (56) implies that k + N > —2 since otherwise k + N would diverge. 
From Eq. (8d) we see that if k > — Vl — 2Ar'^ for some fo, then k will 
increase towards zero and stay positive afterwards. This is, however, 
impossible since then {rN)' = ktN — 2rU'^ < and rN would not converge 
to zero as required. We thus obtain the bounds —2<k< — \/1 — 2Ar'^. 

We then show that W cannot have zeros for r < 1/2. Assume there were 
such a zero at tq. This would imply > 9/4 for r > tq as long as 
\W\ < 1/2. Thus \ W\ = 1/2 for some r = fo + A since otherwise N would 
diverge due to Eq. (11). We then estimate, repeatedly using the Schwarz 
inequality 



and thus N would diverge due to Prop. 4. 

Since the extrema of W are separated by zeros, W is monotonic for suffi- 
ciently large r and therefore has a limit. This limit can only be ±1 since 
otherwise N would diverge due to Eq. (11). 

Next we estimate (using TUW < 0) 



i.e., \TU\ decreases exponentially in r. Eq. (64) then implies that T and 
U both must decreases exponentially since TU cannot change sign. 




(63) 



\TU\-= -2\W\{T^ + 2U^) - k\TU\ < -{y/S\W\ - 2)\TU\ 



(64) 
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Applying Lemma 14 to Eqs. (8d) and (11) we finally obtain k — > — 1 and 
AT — > 1 for r — oo. Thus the solution tends to a regular origin. ■ 



Proposition 6: Suppose a solution of Eqs. (8,9) for a = +1 withN unbounded, 
then the solution has an origin r = at some Rnite t\ with 

i) the S type singularity as described in Eqs. (28,29), if k is unbounded; 

a) the RN type singularity as described in Eqs. (22,23), if k is bounded. 

Proof: In view of Prop. 4 N must be unbounded from below. Integrating 
Eq. (11) shows that A^ — oo for some finite ri, and we assume ri = without 
restriction. In the following we will consider sufficiently small |r| such that 
N(fo) < —1. In order to proceed we distinguish the cases that k is unbounded 
or bounded from below. 

i) For K unbounded from below we first note that k + N is unbounded from 
below and thus k + N ^ — oo for r — due to Eq. (56). 

Prom the equation 

(2A^ - «)•= [N - k){2N -k) + ^(3-3N^ - 2U^ - ST^ - Ar^) 

< (N - k)(2N - k) , (65) 

we then see that if (2 AT — «;)(to) < for some tq, then 2N — k < for all 

T > Tq. 

First assume 2N — k > for all r and therefore k — ^ — oo. Now consider 
the equation 



1 K 



NN 



With A(t) = f^dT'/\N\ = 0(r) we then find (e^K/N)' < 3Ne^ < and 
therefore k/N, which is non-negative, has a finite limit. Lemma 14 then 
implies ndr' = — oo and therefore Ndr' = — oo. This is, however, 
impossible since k/N would then become negative due to Eq. (66). 

Therefore 2A^ — «; < for |r| sufficiently small. Next we want to show 
that K < for |r| sufficiently small. For r ^ this follows from Eq. (8d). 
For r — ^ ri > we first assume that W is bounded. Lemma 13 applied 
to the equation 

(rUy = W(W'^ - 1) + (2N - K)rU , (67) 

then implies that U is bounded and the constraint Eq. (9) implies k < 0. 
Then we assume that W is unbounded and therefore, due to Prop. 4, 
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WU > and \W\ oo. This requires that U is unbounded and the 
equation 

b)= ,3 -(^ + ^)^^ (68) 

then imphes that U jr^ is monotonic and thus IJ diverges. Eq. (8d) then 
again imphes k < 0. 

Now consider 

(tk + rN) = rK{2N - k) + r(l - 2Ar^) . (69) 

The r.h.s. is bounded from below and therefore rn + rN, which cannot be 
positive, has a finite limit. This implies that r — and that vk and rN 
are bounded. Prom the equation 

(r^AT)- = {rK + rN)rN - 2r'^U^ , (70) 

then follows that r^C/^ is integrable and, using the Schwarz inequality, 
that W has a finite limit Wi . Eq. (67) finally implies that W = rU has a 
finite limit Wi. 

Prom Eqs. (13b) we obtain 

m(T) = m(To) - / '-dr' , (71) 

Jrij) ^> 

and therefore lim rm = — VF^^ — (W^^ — 1)^/2. Thus rK and rN have 

finite limits k\ resp. ni satisfying = 2W\ + (VFf — 1)^ and, due to the 
constraint Eq. (9), 2feini = nf + 2Wx - (W^i - 1)^ = 4W^f . 

We want to show that n\ = 0. Assuming ni < 0, we find 2ni — fci < and 
r(r) = niT + O(r^). This implies J^{2N — K)dT' = — oo and Lemma 13 
applied to Eq. (67) yields ki = Wi =0. As a consequence Kdr is 
finite and Eq. (8d) implies that k is bounded from below. This is a 
contradiction, therefore ni = fci = Wi = and Wi = 1, and we may 
assume Wi = +1 without restriction. 

Prom the equation 

{2K-Ny = ^{K+N){2K^N)+6U^+2{1-Ar'^) > -{k+N){2k-N) , (72) 

we see that if (2k — N)(fQ) > for some tq, then 2k — > for all 
To < r < 0. This is, however, impossible because then (2rK — rN)' > 
and 2rK—rN could not converge to zero as required. Therefore 2k— N < 
for all sufficiently small |r|. 

We can now rewrite Eq. (66) as 

/ K\ , , T"^ f K \ Ar^ Ar^ — Ik , ^ 

(-)=(2.-iV) + _(- + l)-_ + _-. (73) 

The first two terms on the r.h.s. are both negative and the last two terms 
are bounded. Therefore 2k — N and T'^/N are integrable and k/N has a 
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limit which must be 1/2 since k + N < 3/r + 0{t) due to Prop. 1. This 
then imphes tN — >■ 2 and Inr/ln |r| — > 2 for r — ^ 0. 

With 

C(r) =J'{'^- ^) tir' = 0(1) , = e^r i , (74) 

such that D = {2N — k)D < 0, we then derive 

^0 

W{t) = 1-J rUdr' , (75a) 

^"■^^^^^ = ^^""^^^""^^ + ^£ " • ^^^^^ 

Starting with W — 1 = 0(1), these equations first imply rU = 0(r) and 
thus W-l = 0(r2), and finally rU = O(r^) and thus W-l = 0(r''). As 
a consequence U and T are bounded, k = 1/r + 0(r), N = 2/t + 0{t), 
and r = -crMr2/2+0(r2) with some aM < as required by Eqs. (28,29). 

ii) For K bounded from below we proceed like in [3], with some modifications 
since we cannot assume k > 1. Let A be a positive constant such that 
K + A > 1. We introduce B = e^'^ and the new dependent variable 
A = W{W'^ - 1) + NU and rewrite Eqs. (8c,e) as 

{BrN)= {k + A)BrN - 2BrU^ , (76a) 
U = X-kU , (76b) 
X = {3W^ -2U^ -1)U . (76c) 

We want to show, that the r.h.s. of these equations stay bounded as r ^ 0. 
In the following e is an arbitrary, suitably small positive number. Since 
{BrNj < we see that and \BrN\~^ are bounded and monotonic and 
therefore their derivatives are absolutely integrable at r = 0. In particular 
we find 

J r^\N\dT <oo, J ^^dr <oo, J -^^-^-^^Jr < oo . (77) 

Next we want to show that W\BrN\~^ is bounded. In view of Prop. 4 this 
is obvious except when > 1 and WW > 0. In that case wc get, using 
the Schwarz inequality, \U\\BN\^'^dT < oo and hence \W\\BrN\~^dT 
is finite for < e < 1. Integration by parts yields 

W , , r W / , 2C/2n ^ , 



\BrN\ 



w r w , 



with some suitably small negative tq. Since the r.h.s. is bounded and both 
terms on the l.h.s. have the same sign W\BrN\~'' is bounded. 



25 



The constraint Eq. (9) implies that {r{N + A))'^ — {W"^ — Vf is bounded 
from above and hence W and with it rN are bounded. As a consequence 

wc get r(0) = 0. 

Next we use the equation 

{r^U) = a^^^^f^+({l + e)N -i^r'U . (79) 

From Lemma 13 we conchide that r'^U is bounded and consequently |C7|" 
is integrable for any n > on the interval (ro,0). Hence A and thus U 
have a finite limit at r = 0, implying that also n is bounded from above. 
We thus get the RN singularity of Eqs. (22,23). ■ 

Since Props. 3, 5, and 6 exhaust all possible eases, we obtain the following: 

Theorem 7 (Classification for a = +1): A solution of Eqs. (8,9) for a = +1 
starting at some regular point tq (which is not a f.p.) can be extended in both 
directions to one of the following singular points: 

i) a regular origin r = at infinite r as in Eqs. (21 ); 

a) a regular horizon of the black hole or cosmological type at some finite ti 
as in Eqs. (35,36); 

Hi) the RN f.p. at infinite r as in Eqs. (40), where ^ for r — > ±oo requires 
W = 0; 

iv) r = at some finite T\ with the S type singularity as in Eqs. (28,29); 
v) r = at some hnite t\ with the RN type singularity as in Eqs. (22,23). 



4.2 Solutions with <t = — 1 

Due to the presence of the cosmological constant the asymptotics for large r is 
deSitter instead of Minkowski and this requires a = —1. We shall again give a 
complete classification of all possible solutions starting at some regular point. 

Proposition 8: Suppose a solution of Eqs. (8,9) for a = —1 with N{t) > 
and unbounded, then r ^ oo for some finite t\ with the deSitter asymptotics 
of Eqs. (37,38). 

Proof: The constraint Eq. (9) implies 2kN > — 1 and therefore 

- 1 lA^I - 1 

\N\ > 1 implies kN > and \k\ > ^j^yp > -^^2 ' ^^^^ 

i.e. K is unbounded from above. Eq. (56) implies that (n + N)/r is bounded 
from above and therefore r — ^ oo. We then introduce N = N/r and R = n/r 
and see from the constraint Eq. (9) that N and R are bounded from below by 
positive constants. This implies that r ^ ti stays finite as r ^ oo. In order to 
be able to use E/r'^ as a Lyapunov function, we introduce a new independent 
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variable p with dp = rdr such that p —>■ oo for r oo, denote the derivative 
with respect to p by a prime, and obtain 



According to Lemma 11.1 of [19] {E/r^)' vanishes on the hmit points of the 
solution, showing that U/r and W/r tend to zero for p — oo. For TV we find 



iV' = l(A-3iV^-l + ^:^) . (82) 



Applying Lemma 14 we obtain N -^A/3 and further from Eq. (9) k 
2VA73. 

Using s = 1/r as independent variable and introducing the new dependent 
variables = iW^^Ueim^) where = s'^W , Ue = s'^'^UN, rrie = s^^m, and 
< e < 1, we can rewrite Eqs. (38) as 

s^W, = eW,-s'-^^^, (83a) 
sj-U, = 2eU, + s^-^[w,{W^-s'^)+s^^'-^^^) , (83b) 

with iV^ = A/3 — + 2s^(^~^)m£. The dynamical system Eqs. (83) has three 
stable modes with real eigenvalues > e. Thus s~'^y^ is bounded near s = for 
any solution vanishing at s = and rj < e (see e.g. [18]). Assuming that is 
bounded and vanishes at s = for some e = eo > 0, which we know to be true 
for eo = 1, we therefore conclude that the same is true for all e > 2eo/3, and 
thus by induction for all e > 0. For e < 1/4 we can finally convert Eqs. (38) into 
integral equations with convergent integrals, i.e. W, UN, and m are bounded 
as required by Eqs. (37). ■ 

Proposition 9: Suppose a solution of Eqs. (8,9) for a = —1 with N and k 
bounded for all r > tq, tiien the solution runs into 

i) theSdSfp. described by Eqs. (39) (with Kg = —1 only possible if W'^ = 1); 

ii) the RN f.p. described by Eqs. (40). 

Proof: The constraint Eq. (9) implies that r, U, and T and thus W are 
bounded. Moreover r has a limit ri > l/\/2A since otherwise k + N would 
diverge due to Eq. (56). Prop. 2 then implies that the solutions runs into 

i) the SdS f.p., as discussed in Section 3. For reasons given there the case 

= — 1 is only possible if W'^ = 1. 

ii) the RN f.p., as discussed in Section 3. ■ 
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Proposition 10: Suppose a solution of Eqs. (8,9) for a = —I with N un- 
bounded and N{fo) < for some tq, then the solution has an origin (r = 0) at 
some finite ri. Depending on the behaviour of N/k as t ^ ti the solution has 

i) an S singularity as described in Eqs. (29), if N/k —>■ 2; 

a) a pseudo-RN singularity as described in Eqs. (25), if N/k — 1/2; 

Hi) unbounded oscillating N, U, and T as described in [16], if N/k has no 
limit. 

Proof: From Eq. (80) we find that k+N is unbounded from below and tiierefore 
diverges to — oo at some finite ri due to Eq. (56). Without restriction we may 
assume ti = 0. Moreover Eq. (80) imphes 

K- -I- AT -I- 1 

K< for k + N<-1, (84) 

o 

and consequently k diverges as well. Lemma 14 applied to Eq. (8d) then yields 
K < 1/t + 0(t). In the following we can assume k < —1, restricting |r| to 
sufficiently small values. 

Next we introduce new dependent variables n = N/k, u = U/k, and t = T/k 
and note that n, u, and t are bounded due to Eqs. (84) and (9). If n has a limit 
ni > then J^^N dr = — oo and r — 0. 

i) First consider the case that n — 2 as r — 0, and consequently u ^ 
and t — > due to the constraint Eq. (9). Prop. 1 then implies tk ^ 1, 
tN 2, and lnr/ln|T| ^ 2 as r — 0. This implies rK 0, rU — > 0, 
and W'^ — > 1 and we may assume W ^ +1 without restriction. 

Defining 

(y -'^)dr' , D = e'^ri , (85) 

such that r'^e'-' is bounded for e > (compare Eq. (74)), we proceed as in 
the proof of Prop. 6 and find that the solution has all properties required 
by Eqs. (28,29). 

ii) Next consider the case that n has a limit ni < 2. We introduce a new 
independent variable p with dp = \K\dT such that p — > oo for r — > and 
denote the derivative with respect to p by a prime. We then obtain again 
Eqs. (30) with z = 1/\k\. 

First assume rii = and therefore u ^ and t ^ 0. Eqs. (30c, c) then 
first imply that has a finite non- vanishing limit and that e'^t and thus 
T are bounded. Eqs. (30b, d) then imply that e^^u and e^^n are bounded, 
i.e. N ^ contrary to the assumption that N is unbounded. 

For < ni < 2 we see from the constraint Eq. (9) that t^ + 2u^ has the 
limit ni(2 — ni) > 0. From Eq. (30b) we see that for any e > and p 
sufficiently large u is monotonic as long as |2u^ — ni\ > e and \u\ > e. 
Since u is bounded it therefore has a limit ui with (2ti^ — ni)ui = 
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and thus t has a hmit ti. Prom Eq. (30d) we then obtain the condition 
2ui{ni + 1) = ni(2 — ni). This only leaves the possibility ni = \ui\ = 
\ti\ = 1/2. Prop. 1 then implies tk — > 2, tN — > 1, and lnr/ln|r| ^ 1 as 
r — 0. Thus r\U\ is integrable and W has a limit Wi, and consequently 
—rN, —rK/2, and \rU\ all have the limit \Wi — 1|. 

Integrating Eq. (8b) we then find W = Wi + 0{t) and therefore ^ 1 
since otherwise f ^ for r ^ 0. Thus the solution has a pseudo-RN 
singularity as described in Eqs. (25). 

iii) Pinally assume that n has no limit. Since n is bounded this implies 
that n has an infinite sequence of maxima and minima where the r.h.s. 
of Eq. (30d) changes sign. Therefore n, u, and t vary between their 
respective maxima and minima. 

If n dp' is finite then the r.h.s. of Eq. (30d) is integrable since < 
n(2 — n) + z'^ due to the constraint Eq. (9), and therefore n has a limit 
contrary to the assumption. Thus N dr' = — n dp' diverges and 
r ^ 0. ■ 

Remark: The proof of Prop. 10 also applies to the case A = corresponding 
to Bartnik-McKinnon type black holes inside their horizon discussed in [16]. 
In fact, this is to be expected, because the cosmological constant becomes less 
relevant as r decreases. 

Since Props. 3 and 8-10 exhaust all possible cases, we obtain the following: 

Theorem 11 (Classification for <t = —1): A solution of Eqs. (8,9) for 
a = —1 starting at some regular point tq (which is not a f.p.) can be extended 
in both directions to one of the following singular points: 

i) r = oo at some finite t\ with the deSitter asymptotics of Eqs. (37,38); 

a) a regular horizon of the black hole or cosmological type at some finite ri 
as in Eqs. (35,36); 

iii) the SdS f.p. at inGnite r as in Eqs. (39), where Kg = Tl for r — ±00 
requires = 1; 

iv) the RN f.p. at infinite r as in Eqs. (40); 

v) r = at some finite ti with the pseudo-RN type singular behaviour of 
Eqs. (25,27); 

vi) r = at some Enite t\ with the S type singular behaviour of Eqs. (28,29); 

vii) r = at some finite ti with unbounded oscillating N, U, and T as de- 
scribed in [16]. 

The cases ii), v), and vi) correspond to the NARN, NAPRN, and NAS solutions 
of [16]. 
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5 4d Space-Times 



In the previous section we have determined aU possible Id radial solutions 
extending between singular points. Taking into account the time independence 
of the solutions this gives us pieces of 2d space-times with a metric as given in 
Eq. (2). Our aim in this section is to determine all inextensible (in the sense of 
geodesic completeness) 4d space-times, which can be constructed from these 2d 
pieces. A standard way for that is to construct the Carter-Penrose diagrams of 
these extensions. Recipes for their construction are given in [23, 24]. As a first 
step, we have to bring the 2-metric to the required standard form 

dsl = Fdt^ - jdr^ , (86) 

by the coordinate change df = re'^dr, where F = ae^^ in terms of our variables. 
Introducing Eddington-Finkelstein type coordinates via du± = dti: F~^df this 
takes the form 

dsl = Fdu+du_ , (87) 

useful for a Kruskal type extension. 

The next step is to express F as a function of f and study its behavior near 
the singular points. This will be done one by one: 

• At r = the solutions are geodesically incomplete. The singular solutions 
have a genuine curvature singularity and do not allow for an extension. 

Although solutions with a regular origin are easily extended taking into 
account the angular dependence, we treat them as inextensible in the 2d 
setting. 

• At a regular horizon the space-time is geodesically incomplete in the space- 
like as well as in the light-like direction and F has a simple zero at r = f/j. 

• At the RN and SdS fixed points the space-time is always inextensible 
(i.e. geodesically complete) in the space-like direction. While for Kg ^ 
the coordinate f tends to infinity for r — ^ ±oo it stays finite for Kg ^ 0. 
Since f plays the role of an afhnc parameter for light-like geodesies, the 
space-time is geodesically complete in the first case and incomplete in the 
second one. For kt — — oo one easily finds f -|- e'^^'^Tg/ i^s and since 

F = (7e2-«(7%(f,-f)2 (88) 

near the f.p. This shows that for ^ the space-time has a degenerate 
horizon. 

There is an exceptional situation of the RNdS f.p. for Hs = 0, which is only 
possible for A = 1/4, where F has a cubic zero at r = \/2. Again one has 
to distinguish the ^ and K ^ for r ^ ibcx), leading to geodesic 

completeness resp. incompleteness at this singular point. For a = +1 
only the RNdS solution with = is possible, but for cr = — 1 there 
are solutions with non-trivial W . They are C°°, but have an essential 
singularity at the degenerate horizon. 



30 



• For solutions with dS asymptotics the space-time is geodesically complete 
at r = oo. 

Before constructing the Carter- Penrose diagrams we present in Figs. 6 and 7 
a graphical representation of all possible compositions of the pieces of 2d space- 
times discussed in the previous section. In principle there is an infinite number 
of such possible compositions, because one can always insert additional pairs of 
regular horizons. We have, however, restricted ourselves to compositions which 
are compatible with the number of parameters of a generic solution. Thus we 
take into account the necessity to suppress the divergent modes at every singular 
point. For example at every horizon we loose one parameter. In this counting 
the singular solutions at r = axe considered as f.p.s without divergent mode 
(compare Sect. 3) as are the solutions with dS asymptotics. For a given A the 
five autonomous Eqs. (8) are reduced to four taking into account the constraint 
Eq. (9). This results in a 3-parameter family of orbits. Hence we can suppress 
maximally four divergent modes varying also A, implying e.g. that we can have 
a maximum of four horizons. This docs, however, not necessarily imply that 
all cases are realized, i.e. can be actually obtained numerically. Nor does it 
mean that solutions with more than four horizons are excluded. It just means 
that we see no way to construct them varying the available parameters at hand. 
Actually we have not even found any solution with four horizons. 

The graphs in Figs. 6 and 7 are more or less self explanatory, but some 
remarks are in order. We ignore the presence resp. position of an equator, which 
would just lead to a much larger number of diagrams without a particular gain 
in information. Likewise we ignore the type of the singularity at a singular 
origin. 

Diagram 2 represents the generic SdS solution, while 3 corresponds to the 
degenerate case. Diagram 5 is the case of SdS with negative mass. Diagram 6 
is the generic RNdS solution, 9 corresponds to the most degenerate RNdS so- 
lution for A = 1/4 with a cubically degenerate horizon. Diagrams 4, 10, and 16 
represent generic solutions. The space-times 10-23 without horizons either have 
spatial sections with 3-sphere topology or have infinite throats with dS resp. 
AdS type asymptotic behaviour. Diagrams 3 and 8 represent the degenerate 
SdS resp. RNdS solutions as well as ones with non-trivial YM field. Diagram 9 
describes besides the most degenerate RNdS solution also solutions joining the 
a = +1 piece of the latter to a non-trivial a = —1 part at the cubically degen- 
erate horizon. The same property is shared by the diagrams 32, 33, and 42. 

Diagrams 15, 22, and 23 describe the degenerate space-time of Eqs. (20a). 
Diagram 41 describes a solution with a degenerate horizon running into the SdS 
f.p. Such a solution could occur as a limit of a 1-parameter family of solutions 
of the type 8 with YM hair. 

The solutions with a regular origin of [9] are given by diagrams 1, 28, 29, 
and 13 in the order as they occur with increasing A. Replacing the regular 
origin by a regular black hole horizon one obtains the diagrams 2, 36, 34, 
and 29. Replacing the origin by two horizons as in the RN black hole one 
obtains the diagrams 6, 44, 43, and 39. The diagrams 19, 20, and 21 are new 
types of regular solutions starting at the RN resp. SdS f.p. and running into 
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Table 1: Parameters of regular solutions 



n 


Ac 


rc 




r* 


Ar 




be 




6* 




br 


1 


0.33049695 


1.729083 


0.33448289 


1.706678 


0.75000000 




0.37534807 


-0.85135 


0.37417940 


-0.85218 


0.25000000 


2 


0.23896851 


1.934703 


0.25004507 


1.812673 


0.36424423 




0.54743193 


+0.45963 


0.53831284 


+0.48411 


0.42959769 


3 


0.23712384 


1.832376 


0.24706870 


1.686850 


0.29321764 




0.57517232 


-0.23626 


0.56568105 


-0.28415 


0.50882906 


4 


0.24027732 


1.746936 


0.24790151 


1.606049 


0.27032753 




0.57699432 


+0.13396 


0.56950069 


+0.18872 


0.54048970 


5 


0.24275653 


1.688873 


0.24861452 


1.556666 


0.26089512 




0.57578215 


-0.08158 


0.56992909 


-0.13598 


0.55401931 


6 


0.24445557 


1.648272 


0.24907014 


1.524611 


0.25638498 




0.57448821 


+0.05216 


0.56981254 


+0.10345 


0.56042575 



the SdS f.p. or to infinity. 

Figs. 8-13 sliow the elementary building blocks of the Carter-Penrose dia- 
grams corresponding to the diagrams of Figs. 6 and 7, which have at least one 
horizon. The complete diagrams are obtained from those by periodic repeti- 
tion, gluing them along suitable common horizons [23, 24]. As usual space-like 
boundaries are represented by horizontal lines, time-like ones are vertical. 

6 Numerical Results 

In this section we present a selection of numerical results, some of which we 
would like to emphasize. We perform a quantitative study of the limit n — ^ oo 
of the 3-sphere solutions. We present a whole new class of solutions with two 
and three regular horizons. We give a complete description of the phase space of 
solutions for a = +1 with a regular origin. An analogous description for a = —1 
exhibits some important aspects. One cannot expect a complete description due 
to the very complicated structure of solutions with a singular origin caused by 
unbounded oscillations [15, 16]. 

6.1 Solutions with a regular origin and a regular horizon 

In Fig. 14 we show parameters of solutions with a regular origin and a regular 
horizon with n = 1 and 2 zeros of W as well as for the limiting solutions (n = oo) 
connecting the origin and the RN f.p. In Table 1 we present parameters of these 
solutions for n < 6 and some special values of A. When continued beyond the 
horizon into the region with a = —1 they are asymptotically deSitter for A up 
to an n-dependent critical value Ac. For A = Ac, indicated by the points Cn in 
Fig. 14 they run into the SdS f.p. while for A > Ac they have an equator and 
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Table 2: Parameters of 3-sphcrc solutions 



n 


A 


h 


''max 


1 


0.7500000000 


0.2500000000 


1.4142136 


2 


0.3642442330 


0.4295976874 


1.4355566 


3 


0.2932176375 


0.5088290567 


1.5213897 


4 


0.2703275292 


0.5404896987 


1.4826529 


5 


0.2608951239 


0.5540193079 


1.4829530 


10 


0.2512790802 


0.5674090979 


1.4396422 


20 


0.2501165095 


0.5688910784 


1.4224344 


30 


0.2500264064 


0.5690006373 


1.4181848 


40 


0.2500089730 


0.5690215671 


1.4165412 


50 


0.2500038394 


0.5690276981 


1.4157401 


GO 


0.25 


0.5690322659 


V2 



end at r = 0. The equator coincides with the horizon for A = A^< (the open 
circles in Fig. 14) and occurs for cr = =Fl if A ^ A^,. The radius of the horizon 
decreases as A increases and vanishes for A = A^, the points Rn, corresponding 
to the regular 3-sphere solutions. 



6.2 Regular 3-sphere solutions 

Table 2 contains the parameters of some regular 3-sphere solutions with up to 
50 zeros of W. As n increases the values A„, 6„ rapidly converge to Aqo = 1/4 
and 6oo ~ 0.569032, i.e. the point i?oo of Fig. 14. The corresponding limiting 
solution connects a regular origin and the RN f.p. with A = 1/4 and asymptotic 
behaviour (compare Sect. 3.8). 

W^^^rp^. (89) 
r-To {t — tqY 

Numerically we find 6 ~ 1.7884 and tq ~ —3.8717 (normalizing r such that 
r — Inr ^ as r ^ 0). Solutions with (Aqo + SA,boo + Sb) for small 6A, 6b 
will miss the RN f.p. but yield a regular 3-sphere provided k, N, and W or 
U simultaneously vanish for some r. Assuming that the phase of W remains 
essentially unchanged this requires that k and vanish at t + 6 {n + l)7r/2 
in order to have n zeros of W. For solutions starting from the origin we first 
have to choose 6b « —1.26 A to suppress the unstable mode (with eigenvalue 
-|-6) near the f.p., where the solution approximately obeys 

k = f-K^ , (90a) 
r = -2N , (90b) 

N = -—-6A, (90c) 

with the consequence that f^/12 + 6Af — N"^ is constant. Moreover Eqs. (90) 
are invariant under the rescaling (dr, f, N, 6A) — > {dr/X, Xk, A^f, X^N, A^JA). 
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Starting from some ti with k = N = and r = — 1 we can thus integrate 

k = f — , (91a) 



+ 46A{f + 1) , (91b) 

towards decreasing r, adjusting ri and (5A such that k and f diverge at r = tq 
with f/K^ ^ 3/4. Numerically we find 5A « 0.81575 and n - tq 3.8288 
and thTis obtain for A„ and f„, the value of f at the equator, the asymptotic 
expressions 

(^{n + 1)^-9- To) x/A„ - 0.25 ^ 3.8288v^0.81575 3.6387 , (92a) 

f 1 

, " ^ ^ -1.1072 , (92b) 

VA„ - 0.25 ^0.81575 

in good agreement with the numerical results shown in Fig. 15 for n < 50. 

6.3 Solutions with two or more regular horizons 

Fig. 16 shows parameters for solutions with two horizons for a = +1 and various 
values of A. In the limit r/j — > for one or both horizons they reproduce 
the solutions with an origin and a horizon and the regular 3-spheres discussed 
above. The dashed-dotted curve delimits the region of existence; on it the value 
of Wh vanishes and the solutions bifurcate with Reissner-Nordstr0m-deSitter 
solutions. The dotted curves correspond to the critical A values, giving solutions 
running into the Schwarzschild f.p. in the region a = —1. The solutions in 
the region outside the dotted curves are asymptotically deSitter. The dashed 
curves correspond to solutions with coinciding equator and horizon. Solutions 
in the region inside the dashed curves have an equator between their horizons 
in contrast to those outside of it. The RN-deSitter and the critical solutions 
have no equator and thus the dashed-dotted and the dotted curves must lie 
outside the dashed ones. Since on the diagonal the two horizons are equal, these 
solutions have an equator and thus the diagonal lies inside the dashed curves. 
Thus the dashed curves must go through the point P on the intersection of the 
diagonal with the boundary curve, where the two horizons merge. Since the 
solutions on the boundary curve are asymptotically deSitter, the dotted curve 
must also go through P. At this point the solution can be obtained exactly 
through a suitable rescaling. Putting 

W = eW,U = eU,N = e^N and r = R{l + e^f) , (93) 

where e is a small parameter we obtain in leading order from the constraint 
Eq. (9) A = (i?2 - l)/r-4 and from Eqs. (8) 



W = RU (94a) 
W 

U = -—-kU (94b) 
R 

k=l- 2Ai?2 - «;2 (94c) 
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Table 3: Some solutions with three regular horizons and two zeros of W 



A 


Wo 


ri Wi 


r2 W2 


rs W3 


Woo 


0.0045 


1.3655 


0.0948 1.3582 


1.0256 0.1000 


24.6488 0.0862 


0.1637 


0.0643 


1.4072 


0.1207 1.3964 


1.0459 0.0800 


5.6355 0.0263 


0.0638 


0.1187 


1.4388 


0.1417 1.4249 


1.0918 0.0600 


3.7151 0.0197 


0.0610 


0.1655 


1.4737 


0.1664 1.4556 


1.1502 0.0400 


2.8234 0.0152 


0.0647 


0.2061 


1.5141 


0.1968 1.4904 


1.2351 0.0200 


2.2350 0.0098 


0.0716 


0.2247 


1.5378 


0.2155 1.5103 


1.3051 0.0100 


1.9718 0.0059 


0.0762 


0.2337 


1.5511 


0.2263 1.5214 


1.3623 0.0050 


1.8255 0.0034 


0.0789 


0.2411 


1.5630 


0.2361 1.5312 


1.4544 0.0010 


1.6627 0.0008 


0.0813 



Solving Eq. (94c) we get k = a cot (ar) as in Eq. (20c) with = 2Ai?^ — 1 and 
plugging the result into Eqs. (94a,b) we arrive at the second order equation 

W + a cot(ar)T^ + W = . (95) 
Introducing x = cos (ar) we find the Legendre equation 



— {{l-x^)—]+^W = 0. (96) 



The acceptable solutions are those regular at x = ±1 given by the Legendre 
polynomials, obtained for a""^ = n{n + 1) = 2, 6, 12, ... yielding 

, 2 12 24 , , 

B^ = 5^ = 4, -, -, .... (97a) 

l-o^^^^ 143 

4 16' 144' 576' ^ ' 

The cases n = 1 and 2 are shown in Fig. 16. 

Combining numerical results for a = +1 and u = — 1 we have found a one 
parameter family of solutions with three regular horizons and two zeros of W. 
In Table 3 we present the parameters of some of these solutions. Starting from 
r = with an RN type singularity and W = Wq, they have two regular horizons 
at r = ri_2 with W = Wi^2, two zeros of W, a third horizon at with W3, and 
extend to r = 00 with dS asymptotics and W = Woo- There is a similar such 
family with three zeros of W, but we found no such solution with only one zero. 



6.4 Description of the phase space 

Fig. 17 schematically shows the phase space for solutions starting from a regular 
origin, i.e. with a = +1, dividing the (A, b) plane into various regions according 
to their type and including results for A = from [3] resp. for A < from [12]. 
All the cases listed in the Classification Thm. 7 are present in the half plane 
A > of Fig. 17. The regions AdS„ with A < and S„ correspond to solutions 
with n zeros of W and AdS asymptotics or a singular origin respectively. At 
the moment this picture is mostly based on numerical results, but proving the 
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existence of such a picture will certainly be an essential part in an existence 
proof for the various solution types. 

In order to simplify the description of the phase space we deform and com- 
pactify Fig. 17 into Fig. 18 as follows: First we identify the points (A, ±00) and 
add a dashed-dotted line along b = ±00 separating the two previously discon- 
nected parts of the region Sq- Next we add a point Rq = (+00, 0) corresponding 
to the limit r/j ^ of deSitter space and compactify (+00, b) to this point. The 
points with constant (finite) A are then on concentric curves around Rq. We 
then add the points P_i, and R-i as the limits A —>■ —00 of (A, — y^|A|/4), 
(A, +Y^|A|/4), and (A, ±00). And we add points P„ for n = 1, ... all identical to 
Poo- Finally the boundary of Fig. 18 corresponding to (—00,6) is compactified 
to one point 6_i = P_i = P_i. 

The open region AdS = U„>oAdS„ is the complement of the closure of the 
open region S = U„>oS„. The points 6„ at (0,6„) represent the asymptotically 
flat BK solutions (augmented by Minkowski space for n = 0), with their limit 
point (compare [3]). The points P„ represent the regular 3-sphere solutions, 
with their limit point Poo at A = 1/4. AdS„ and AdS„+i are separated by the 
dotted curves from 6„ to P„ where W^o = 0; S„ and S„+i are separated by 
the dotted curves from P„ to P„ where Wq = and the solid curves from 
bn to Rn with a regular horizon. The dashed-dotted curves from Pn-i to P„ 
with I Wo I = 1 and an S type singularity at the origin divide into parts with 
\Wo\ > 1 and |VFo| < 1, both with a RN type singularity. Similar dashed-dotted 
curves from bn-i to bn where \ Woo\ = 1 separate AdS^ into parts with |VFoo| > 1 
and I Woo I < 1- The regions AdS„ and S„ are separated by the solid curves from 
bn-i to bn with SAdS asymptotics where Wx) resp. Wq diverge (compare [12]). 
Finally the dashed curve from Poo through b^o, Poo with A = —15/4, and Pq 
to P_i represents solutions with the RN f.p. asymptotics. They oscillate for 
A > —15/4 and therefore all regions AdS„ and S„ for n S> 1 are in the vicinity 
of the curve from Poo to Poo- For A < —15/4 the RN solutions have a finite 
number of zeros; numerically we found n = 1 from Poo to Pq with A ~ —5.0646 
and n = from Pq to P_i. The curves of constant Woo in the region AdS„ 
all connect bn-i and 6„ whereas the curves of constant Wq in the region all 
connect Pn-i and P„, approaching the solid boundary as Woo resp. Wq diverges 
or the dotted and dashed one as Woo resp. Wq vanishes. 

Fig. 19 shows the phase space of solutions for a = —1 starting from a regular 
horizon with at rh = 0.1, again with all the cases listed in the Classification 
Thm. 11. The open region dS = U„>odSn of the (A, W^) plane corresponds 
to solutions with dS asymptotics; this region is divided into regions dS„ with 
n zeros of W, separated by the dashed-dotted curves where H^oo = 0. The 
boundary of dS represents solutions ending in a f.p. with finite r: the SdS f.p. 
with W = (—1)" at the solid curves bounding dS„ or the RN f.p. with W = 
at the points Qi connecting these curves. 

The complement of the closure of dS represents different types of solutions. 
They can end at a second regular horizon (dashed curves) or a singular origin 
with either an S type singularity (dotted curves) , a pseudo-RN type singularity 
(cross at (A,W/i) ^ (0.09,0.5)), or unbounded oscillations. In Fig. 19 only a 
selection of such curves is shown. 
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7 Appendix 



First we extend Prop. 1 of [3] by a uniqueness statement: 

Proposition 12: Suppose a system of differential equations for m+n functions 
u = (ui, . . .,Um) and v = {vi, . . .,Vn), 

= f^fiit, u, v) , t^ = -XiVi + t^^giit, u, v) , (98) 

with constants Aj with Re(Aj) > and integers f-iijUi > 1 and let C be an open 
subset of ET^ such that the functions / and g are analytic in a neighbourhood 
of t = 0, u = c, V = for all c e C. Then there exists a unique m-parameter 
family of solutions of the system (12) with boundary conditions {u,v) = (c, 0) 
at t = 0, defined for c^C and \t\ < tQ{c) with some to{c) > 0. These solutions 
satisfy 

Ui{t) = ci + Oit^"') , Vi{t) = Oit"^) , (99) 
and are analytic in t and c. 

Proof: The proof of existence and analyticity in [3] for real Aj > applies 
without any change to complex values with positive real part. 

In order to prove uniqueness we assume (without loss of generality) t > 0, 
introduce r = — Ini as new independent variable as well as s = y/t as additional 
dependent variable, and replace Eqs. (98) by the autonomous system 

s = ~, (100a) 

Ui='^-s''^'i'^fi{s^c + su,v) , (100b) 

Vi = XiVi - s'^"^ gi{s'^ , c + su, v) , (100c) 

where Ui = {ui — Ci)/s vanishes as r — > oo, i.e. s — due to Eqs. (99). This 
system has the hyperbolic f.p. s = n = v = with one stable and m+n unstable 
modes. Thus there exists a unique 1-dimensional stable manifold with (n, v) 
analytic in s. Since Eqs. (100) are invariant under the substitution [s,u,v) 
(— s, —u, v) this implies that (u, v) is analytic in t. ■ 

Next we prove some simple Lemmata used in the proof of the classification 
theorems. 

Lemma 13: Consider a solution y of the linear differential equation y = a + by 
in some interval tq < t < ti with \a\ integrable. If 

c(r',r) = J\{T")dT" , (101) 

is bounded from above for tq < t' < t < Ti then y is bounded; if c{t' , r) has a 
limit as T ^ Ti then yir) has a limit; if c{t' , ti) = — oo then y{Ti) = 0. 

Proof: All properties are implied by the explicit form 

y(r) = y{TQ)e<^'^^^ + f a{T')e<^' '^Ut' . (102) 

J TO 
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Lemma 14: Consider solutions of the Riccati equation y = a — y^. 

i) Suppose \a\ is integrable on some interval tq < t < ti, then the solution 
which is unbounded at ti behaves as y = 1/{t — ti) + 0{1) as r — ^ ri. If 
a is bounded then y = l/(r — ti) + 0{t — ri). 

a) Suppose a is bounded from below and y < for tq < t < ti and y diverges 
as r — Ti, then y < l/(r — ri) + 0(r — ri). 

Hi) Suppose a has a limit ai for r — > oo and y is bounded, then ai must be 
non-negative and y{T) — > for r — ^ oo. 

Proof: 

i) Without restriction we assume ri = 0, put y = 1/t + rj, and get by 
integration 




(103) 



which we solve by iteration starting with r)o = 0. Putting 

^ = Iroi^' MMr')\dT' and ||/|| = sup{|/(r)|, tq < r < 0} we can 
estimate ^ ^(l~l'^ol/3) ""^ < 1 choosing |ro| small enough. From this 
we get ||??n+i ~%|| ^ (2|'^o|/3)"||?7i|| and thus convergence for |ro| < 3/2. 
For \a\ < B and TqB < 1 we can estimate ||?7n/T|| < B/2. 

ii) Let A be a constant such that a > A. We assume again ri = 0, put 
r) = l/(ry), and get by integration 

r){T) = 1 + y —av'^ dr' <1 + A J —rf' dr' . (104) 

For A > we obtain immediately r/ < 1 and thus y < 1/r. For A < Q 
we choose |ro| sufficiently small such that I^Itq < 1/2 and can estimate 
77 < 1 + \ A\t'^ and consequently y < 1/r — \A\t. 

iii) Obviously ai must be non-negative, otherwise y would be unbounded from 
below. For any e > we can find such that |a — ai| < e/2 for r > r^. 
If y(Tc) < — "v/ai + e, then y will be unbounded from below. Now suppose 
that y{Te) > \/ai + e, then we can find r^' > such that y decreases for 
Te < T < Tg and y < + e for r > r^'. If ai = this implies y ^ 
for r — 00. For e < ai 7^ now suppose that y^{Te) < al — e, then we 
can find r^' > such that y increases for Tg < r < r^' and y > — e 
for T > r^. This implies y ^ ^/o? for r ^ 00. On the other hand, if 
— ^/a^+^ < y{T^) < — V^i ~ £ for any e then y ^ — ^ 



Lemma 15: Consider a solution of the Riccati equation y = a + hy — y"^. If a 
is bounded from above and b is bounded for t > tq, then y is either bounded 
for all T >To or diverges to —00 for some Rnite n > tq. 
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Proof: Let A, B be positive constants such that a < A and \b\ < B. We can 
estimate y < for |y| > C = ^/2A + 2B and therefore y is bounded from above. 
Furthermore y monotonically decreases and > 1/2 for y < —C, and thus 

y — — GO for some finite ri. ■ 
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Figure 1: Solutions near a singular origin with |Wo| ~ 1 
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Figure 2: Four trajectories of Eqs. (32) emerging from the focal point 77, () 
(1/2, 1, 1/2) and rapidly approaching the bounding triangle 
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Figure 5: RN f.p. with A = 1/4: Tangential flow after the blow-up 
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Figure 7: Id graphs with horizons 
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Figure 8: Carter-Penrose diagrams: one regular or degenerate horizon 




Figure 9: Carter-Penrose diagrams: two regular horizons 




Figure 10: Carter-Penrose diagrams: one degenerate and one regular horizon 
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Figure 11: Carter- Penrose diagrams: three regular horizons 
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Figure 13: Carter-Penrose diagrams: four regular horizons 
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Figure 14: Regular solutions with n = 1, 2, and oo zeros of W 
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Figure 15: Asymptotic behaviour of the parameters of regular 3-sphere solutions 




Figure 16: r/^ vs. r/^j for solutions with n = 1 and 2 zeros of W 
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Figure 18: Compactified phase space of solutions for cr = +1 with a regular 
origin 
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